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Abstract. An overview of recent HERA results on inclusive production of D f± mesons in deep inelastic 
scattering is given. 
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1 Introduction 

The charm mass, m c , is larger than the QCD dimensional 
scale Aqcdi therefore, perturbative QCD is applicable at 
the scale m c . In deep inelastic scattering (DIS), another 
scale is the squared four- momentum transfer, Q 2 , carried 
by the exchange photon. The conventional QCD interpre- 
tation for Q 2 ~ m 2 is that charm is determined solely 
by the gluon density, i.e. charm quarks are generated dy- 
namically through the boson-gluon fusion (BGF) process. 
Such approach is called the fixed-flavour-number scheme 
(FFNS). For a sufficiently high Q 2 , this description may 
not be adequate, thus the interplay between two indepen- 
dent scales, m c and Q, embodies interesting QCD physics. 

Since the charm quarks are copiously produced via 
gluon splitting, charm mesons can be used for testing dif- 
ferent resummation techniques. The DGLAP resumption 
is often used together with the BGF process calculated up 
to next-to-leading-order (NLO) QCD (the HVQDIS pro- 
gram I ). The BGF mechanism at leading-order QCD is 
implemented in parton-shower Monte Carlo (MC) mod- 
els (AROMA g] and RAPGAP [3]). Another description 
of the charm production is based on the CCFM evolu- 
tion 0] as implemented in the CASCADE model 5 . This 
also starts from the BGF process but convoluted with the 
unintegrated gluon density. 

The ZEUS analysis discussed in this overview was per- 
formed with the data (82 pb" 1 ) taken from 1998 to 2000, 
when electrons or positrons with energy E e = 27.6 GcV 
were collided with protons of energy E p — 920 GeV. The 
HI Collaboration uses 1996-1997 data, when HERA op- 
erated with E p = 820 GeV. In this paper, only high- 
statistics results based on reconstructed D*^ mesons are 
discussed (other D mesons are discussed in |H]). The mesons 
were identified using the decay channel D* + — * D°ir + with 
the subsequent decay D° — > K ~n + and corresponding an- 
tiparticle decay. 

a For the HI and ZEUS Collaborations 
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Fig. 1. Differential D* ± cross section as a function of Q 2 for 
0.02 < y < 0.7, 1.5 < pt(-D* ± ) < 15 GeV and \rj(D* ± )\ < 1.5, 
compared to the NLO QCD calculation of HVQDIS using the 
Peterson fragmentation with e = 0.035 and the hadronisation 
fraction f(c — » D*) = 0.235 (solid line). The renormalisation 
and factorisation scales were set to yj Q 2 /4 + m 2 . Also shown 
are the NLO QCD predictions based on the CTEQ5F3 PDF 
(dashed-dotted line) and an alternative hadronisation scheme 
(dotted line). The largest theoretical uncertainties (filled aria) 
are due to the charm mass (±0.15 GeV) and the scale varia- 
tions (0.5l4 H , 2fi R ). 



2 Inclusive cross sections 

Figure 1 shows the differential D*^ cross section as a func- 
tion of Q 2 [7j . The data falls by about four orders of mag- 
nitude in the measured region. The NLO QCD calcula- 
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Fig. 2. Differential D"" ± cross sections as a function of T]{D* ) 
compared to the NLO QCD calculation of HVQDIS. The cross 



section was measured for Q 2 > 1.5 GeV^ 
Fig. 1. 



Other details as for 



tions based on the FFNS give a good agreement with the 
measurements up to Q 2 = 1000 GeV 2 . Predictions using 
an alternate PDF, CTEQ5F3, and the Lund-string hadro- 
nisation extracted from AROMA, instead of the Peterson 
fragmentation, are shown separately. 

The pseudorapidity distribution of D*^ mesons is known 
to be particularly sensitive to the underlying parton dy- 
namics at small Bjorken x jHJ. Figure 2 shows that the 
NLO calculation based on the ZEUS NLO fit [§], together 
with the Lund string fragmentation from AROMA, gives 
the best description of the ^(D*^) cross section (and also 
better than the prediction using GRV98-HO, not shown). 
The NLO predictions with the CTEQ5F3, and with the 
Peterson fragmentation, have distributions which are less 
shifted in the forward region. 

The HI Collaboration also observes differences between 
data and NLO QCD when the Peterson fragmentation and 
the CTEQ5F3 PDF is used for the NLO calculations US], 
Fig. 3. The agreement with the data is better when the 
CASCADE model based on the CCFM evolution is used. 
Since both HVQDIS and CASCADE use the Peterson 
fragmentation, the difference between these models can 
be attributed to the parton dynamics at low x. 

In contrast, ZEUS uses the MC models with the Lund 
string fragmentation. In Fig. 4, AROMA and CASCADE 
are compared with the data. The CASCADE overesti- 
mates the data, while AROMA is slightly below. The theo- 
retical uncertainties, which are expected to be larger than 
for the NLO calculations, were not estimated. At present, 
data are not precise enough to distinguish between the 
shapes of the 7](D* ± ) distributions. 

The resolved processes, in which the photon displays 
a hadronic structure, can contribute to the D*^ cross sec- 
tion at low Q 2 . The resolved contribution, as implemented 
in the RAPGAP model, enhances the D*^ cross section 
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Fig. 3. Differential D cross section for Q > 2 GeV and 
0.05 < y < 0.7 as a function of r](D* ± ) compared with the 
CASCADE and HVQDIS programs. The kinematic range for 
the Pt(D* ± ) as for Fig.l. Both theoretical predictions are 
based on the Peterson fragmentation with the parameters as 
discussed in [10]. 



in the rear direction JSj- Therefore, it is unlikely that the 
observed differences between the NLO QCD and the data 
can be explained by the resolved processes. 

Generally, the description of the total charm cross sec- 
tion by MC models is not perfect. For more detailed tests, 
the total inclusive D* cross section was calculated together 
with the D* cross section with associated dijets JS]j as 
shown in Fig. 5. The CASCADE model is systematically 
above the data, while RAPGAP (LO BGF) is below. 

Figure 4 shows the ratio of e~p and e + p £>* ± cross 
sections as a function of rj(D* ± ) for Q 2 > 1.5 GeV 2 . This 
ratio exhibits a trend to increase with increasing Q 2 and 
x . For Q 2 > 40 GeV 2 , there is ~ 3a difference between 
e~p and e + p D*^ cross sections. According to the Stan- 
dard Model, charm cross sections do not depend on the 
charge of the lepton in ep interactions. This difference is 
assumed to be a statistical fluctuation, and two sets of 
the data were combined. More e~p data from HERA II 
is necessary to show whether the difference with the e + p 
data is indeed a statistical fluctuation. 



3 Extrapolation results 

A popular way to look at double differential charm cross 
sections as functions of Q 2 and x is to reconstruct F| c . In- 
tegrated cross sections in Q 2 and y kinematic bins were ex- 
trapolated to the full phase space using HVQDIS based on 
the Peterson fragmentation function. Usually, several un- 
certainties in the extrapolation are considered. The largest 
uncertainties are those associated with the AROMA model 
for fragmentation and charm-mass variations (±0.15 GeV). 

Figure 6 shows the F^as a function of x at Q 2 values 
between 2 and 500 GeV 2 . The data rise with increasing 
Q 2 ; the rise becoming steeper at lower x. Comparisons 
of the ZEUS measurements with previous results from 
HI ^T] show good agreement. The data are well described 
by the NLO prediction based on the ZEUS NLO QCD fit. 
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Fig. 4. Differential D* ± cross section as a function of r)(D*) 
compared with the AROMA (dashed line) and CASCADE 
(solid line) MC programs. Other details as for Fig.l. 



4 Conclusions 

The production of D*^ mesons has been measured in DIS 
at HERA over a wide kinematic range in Q 2 , from Q 2 = 
1.5 GeV 2 to Q 2 = 1000 GeV 2 . The NLO QCD based on 
the FFNS gives a good agreement with the data up to the 
highest Q 2 range measured. 

At present, no conclusive statement can be made on 
the applicability of the CCFM evolution, since other ef- 
fects related to charm fragmentation and the gluon den- 
sity inside the proton, are shown to affect the D*^ cross 
sections. 

More data from HERA II is necessary to increase the 
precision of the measurements and to extend the kinematic 
range of reconstructed charm mesons. This should allow 
to understand the applicability of the FFNS at high Q 2 
and the CCFM parton evolution at low well as to 

understand whether there is a difference between e~p and 



p D** data. 
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Fig. 5. The total inclusive D* ± cross section versus the D*^ 
cross section with associated dijet compared to MC models. 
The cross section was measured for Q 2 > 2 GeV 2 , 0.05 < y < 
0.7 1.5 < pr(D* ± ) < 15 GeV and \n(D r± )\ < 1.5. For the D** 
cross section with associated dijets, Eti > 4 GeV, Et2 > 3 
GeV and —1 < r/Qets) < 2.5 were applied. 
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Fig. 6. The measured F^ at Q values between 2 and 
500 GeV 2 as a function of x. The curves correspond to the 
result of the ZEUS NLO QCD fit, where the lower and upper 
curves show the uncertainty on the parton distributions for the 
ZEUS NLO fit. 
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